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We analyze the properties of the ACOT scheme for heavy quark production and make use of 
the MS massless results at NNLO and N 3 LO for the structure functions F2 and Fl in neutral 
current deep-inelastic scattering to estimate the higher order corrections. For this purpose we 
decouple the heavy quark mass entering the phase space from the one entering the dynamics of 
the short distance cross section. We show numerically that the phase space mass is generally 
more important. Therefore, the dominant heavy quark mass effects at higher orders can be taken 
into account using the massless Wilson coefficients together with an appropriate slow-rescaling 
prescription implementing the phase space constraints. Combining the exact ACOT scheme at 
NLO with these expressions should provide a good approximation to the missing full calculation in 
the ACOT scheme at NNLO and N 3 LO. 
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I. INTRODUCTION 



A. Motivation 



H1 and ZEUS 



The production of heavy quarks in high energy pro- 
cesses has become an increasingly important subject of 
study both theoretically and experimentally. The the- 
ory of heavy quark production in perturbative Quan- 
tum Chromodynamics (pQCD) is more challenging than 
that of light parton (jet) production because of the new 
physics issues brought about by the additional heavy 
quark mass scale. The correct theory must properly take 
into account the changing role of the heavy quark over 
the full kinematic range of the relevant process from the 
threshold region (where the quark behaves like a typical 
"heavy particle") to the asymptotic region (where the 
same quark behaves effectively like a parton, similar to 
the well known light quarks {u, d, s}). 

With the ever-increasing precision of experimental 
data and the progression of theoretical calculations and 
parton distribution function (PDF) evolution to next-to- 
next-to-leading order (NNLO) of QCD there is a clear 
need to formulate and also implement the heavy quark 
schemes at this order and beyond. The most important 
case is arguably the heavy quark treatment in inclu- 
sive deep-inelastic scattering (DIS) since the very pre- 
cise HERA data for DIS structure functions and cross 
sections form the backbone of any modern global anal- 
ysis of PDFs. Here, the heavy quarks contribute up to 
30% or 40% to the structure functions at small momen- 
tum fractions x. Extending the heavy quark schemes to 
higher orders is therefore necessary for extracting precise 
PDFs and hence for precise predictions of observables at 
the LHC. However, we would like to also stress the theo- 
retical importance of having a general pQCD framework 
including heavy quarks which is valid to all orders in per- 
turbation theory over a wide range of hard energy scales 
and which is also applicable to other observables than 
inclusive DIS in a straightforward manner. 

An example, where higher order corrections are par- 
ticularly important is the structure function Fl in DIS. 
The leading order (0{oP s )) contribution to this structure 
function vanishes for massless quarks due to helicity con- 
servation (Callan-Gross relation). This has several con- 
sequences: 

• Fl is useful for constraining the gluon PDF via the 
dominant subprocess 7*9 — ¥ qq. 



• The heavy quark mass effects of order 0(^ 
relatively more pronounced^ 



arc 



1 Similar considerations also hold for target mass corrections 
(TMC) and higher twist terms. We focus here mainly on the 
kinematic region x < 0.1 where TMC are small [l|j. An inclu- 
sion of higher twist terms is beyond the scope of this study. 
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Figure 1: Fl vs. Q from combined HERA-I inclusive 

deep inelastic cross sections measured by the HI and 

ZEUS collaborations. Figure taken from Ref. Q. 



• Since the first non-vanishing contribution to Fl 
is next-to-leading order (up to mass effects), the 
NNLO and N 3 LO corrections are more important 
than for F2. 

In Fig. [1] we show a comparison of different theoretical 
calculations of Fl with preliminary HERA data [2j . As 
can be seen, in particular at small Q 2 (i.e. small x), there 
are considerable differences between the predictionso 

The purpose of this paper is to calculate the leading 
twist neutral current DIS structure functions F^ and Fl 
in the ACOT factorization scheme up to order 0(a 3 s ) 
(N 3 LO) and to estimate the error due to approximating 
the heavy quark mass terms 0(a 2 s x : g 3 -) and 0(a 3 s x , -£ T ) 
in the higher order corrections. The results of this study 
form the basis for using the ACOT scheme in NNLO 
global analyses and for future comparisons with precision 
data for DIS structure functions. 



B. Outline of Paper 

The rest of this paper is organized as follows. In Sec. 
[II] we review theoretical approaches to include heavy fla- 
vors in QCD calculations. Particular emphasis is put 
on the ACOT scheme which is the minimal extension 
of the MS scheme in the sense that the observables in 
the ACOT scheme reduce to the ones in the MS scheme 
in the limit m — >• without any finite renormalizations. 
In this discussion we explicitly distinguish between the 
heavy quark/heavy meson mass entering the final state 
phase space which we will call "phase space mass" and 
the heavy quark mass entering the dynamics of the short 



2 An updated analysis of the HI measurements extending down 
to even lower Q 2 values has been published in Ref. Q, and a 
combined analysis with ZEUS is in progress. 
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Figure 2: Characteristic Feynman graphs which con- 
tribute to DIS heavy quark production in the ACOT 
scheme: a) the LO 0(a s ) quark-boson scattering QV — > 
Q, b) the NLO 0(a 1 s ) gluon-boson scattering gV — > QQ, 
and c) the corresponding subtraction term (SUB) (g —$■ 
QQ)®{Q^gQ). 



distance cross section denoted "dynamic mass." We show 
numerically using the exact ACOT scheme at 0(as) 
(NLO) that the effects of the phase space mass are more 
important than the ones due to the dynamic mass. We 
use this observation to construct in Sec. Mil the NC DIS 
structure functions in the ACOT scheme up to C(a|). 
The corresponding numerical results are presented in 
Sec. IIV1 Finally, in Sec. |V| we summarize the main re- 
sults. 



II. REVIEW OF THEORETICAL METHODS 

We review theoretical methods which have been ad- 
vanced to improve existing QCD calculations of heavy 
quark production, and the impact on recent experimen- 
tal results. 



A. ACOT Scheme 

The ACOT renormalization scheme |4| provides a 
mechanism to incorporate the heavy quark mass into the 
theoretical calculation of heavy quark production both 
kinematically and dynamically. In 1998 Collins [5J ex- 
tended the factorization theorem to address the case of 
heavy quarks; this work provided the theoretical foun- 
dation that allows us to reliably compute heavy quark 
processes throughout the full kinematic realm. 

Figure [2] displays characteristic Feynman graphs for 
the first two orders of DIS heavy quark production. If we 
consider the DIS production of heavy quarks at 0(a s ) 
this involves the LO QV —t Q process and the NLO 
gV -> QQ process^ 

The key ingredient provided by the ACOT scheme is 
the subtraction term (SUB) which removes the "double 



3 At NLO, there are corresponding quark-initiated terms; for sim- 
plicity we do not display them here, but they are fully contained 
in our calculations 



counting" arising from the regions of phase space where 
the LO and NLO contributions overlap. Specifically, at 
NLO order, we can express the total result as a sum of 

CfTOT = &LO + WnLO — VSUb} (1) 

where the subtraction term for the gluon-initiated pro- 
cesses is 



&SUB = fg® Pg^Q <8> &QV^Q ■ 



(2) 



a sub represents a gluon emitted from a proton (f g ) 
which undergoes a collinear splitting to a heavy quark 
{Pg^o) convoluted with the LO quark-boson scattering 



(Jqv^q- Here, 



P g ^ Q (x,n) = §£ ln(fi 2 /m 2 



P„ 



}(x) 



where P g ^Q(x) is the usual MS splitting kernel, m is 
the quark mass and fi is the renormalization scalqfl which 
we typically choose to be fi = Q. 

An important feature of the ACOT scheme is that it 
reduces to the appropriate limit both as m — > and 
in — > oo as we illustrate below. 



1. Fixed- Flavor- Number- Scheme (FFNS) Limit 

Specifically, in the limit where the quark Q is rela- 
tively heavy compared to the characteristic energy scale 
(/i < to), we find ot,o ~ a sub such that (Jtot ~ °nlo- 
In this limit, the ACOT result naturally reduces to 
the Fixed-Flavor-Number-Scheme (FFNS) result. In the 
FFNS, the heavy quark is treated as being extrinsic to 
the hadron, and there is no corresponding heavy quark 
PDF (/q ~ 0); thus a LO ~ 0. We also have a sub ~ 
because this is proportional to ln(/i 2 /to 2 ). Thus, when 
the quark Q is heavy relative to the characteristic energy 
scale [i, the ACOT result reduces to ijtot ~ <?nlo- 



2. Zero-Mass Variable-Flavor-Number-Scheme 
(ZM-VFNS) Limit 

Conversely, in the limit where the quark Q is rela- 
tively light compared to the characteristic energy scale 
(/i > to), we find that ouo yields the dominant part 
of the result, and the "formal" NLO 0(as) contribution 
{cnlo — csub} is an 0{as) correction. 

In the limit m//x — > 0, the ACOT result will reduce 
to the MS Zero-Mass Variable-Flavor-Number-Scheme 
(ZM-VFNS) limit exactly without any finite renormal- 
izations. In this limit, the quark mass to no longer plays 
any dynamical role and purely serves as a regulator. The 
cnlo term diverges due to the internal exchange of the 
quark Q, and this singularity will be canceled by ijsub- 



4 In this subsection we will distinguish /i and Q; in the following, 
we will set /j, = Q and display the results as a function of Q. 
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Figure 3: Comparison of F${x, Q) (scaled by 10 4 ) vs. the 
quark mass m in GeV for fixed x = 0.1 and Q = 10 GeV. 
The red dots are the full ACOT result, and the blue 
line is the massless MS result. The logarithmic plot 
demonstrates this result holds precisely in the m — > 
limit. 



3. ACOT as a minimal extension of MS 

We illustrate the versatile role of the quark mass in 
Fig. Oa where we display the ~MS ZM-VFNS and the 
ACOT result as a function of the quark mass m. 

We observe that when m is within a decade or two of 
/i that the quark mass plays a dynamic role; however, 
for m <C n, the quark mass purely serves as a regulator 
and the specific value is not important. Operationally, it 
means we can obtain the MS ZM-VFNS result either by 
i) computing the terms using dimensional regularization 
and setting the regulator to zero, or ii) by computing the 
terms using the quark mass as the regulator and then set- 
ting this to zeroO To demonstrate this point explicitly, 
in Fig. H-b we again display the WS ZM-VFNS and the 
ACOT results but this time with a logarithmic scale to 
highlight the small m region. We clearly see that ACOT 
reduces the WS ZM-VFNS exactly in this limit without 
any additional finite renormalization contributions^ 



5 If we were to compute this process in the MS scheme, the 
In (m 2 /Q 2 ) in the SUB term would simply be replaced by a 
1/e pole which would cancel the corresponding singularity in 
the NLO contribution. 

6 It is possible to define other massive schemes that could include 
additional matching parameters or extra observable— dependent 
contributions. For example, the calculation of F^ in the original 
RT scheme [7| included extra higher-order contributions that do 
not vanish as Q/m —¥ oo. 
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Figure 4: Comparison of the DGLAP evolved charm 

PDF f c (x,fi) with the perturbatively computed single 

splitting (SUB) f c (x, //) = f g (x, fi) <8> P g ^ c vs. fj, in GeV 

for two representative values of x. 



The ACOT scheme is minimal in the sense that the 
construction of the massive short distance cross sections 
does not need any observable-dependent extra contribu- 
tions or any regulators to smooth the transition between 
the high and low scale regions. The ACOT prescription 
is to just calculate the massive partonic cross sections 
and perform the factorization using the quark mass as 
regulator. 

It is in this sense that we claim the ACOT scheme 
is the minimal massive extension of the MS ZM-VFNS. 
In the limit m/fi — > it reduces exactly to the MS 
ZM-VFNS, in the limit m/fi > 1 the heavy quark de- 
couples from the PDFs and we obtain exactly the FFNS 
for m//i ^> 1 and no finite renormalizations or additional 
parameters are needed. 



4- When do we need Heavy Quark PDFs 

The novel ingredient in the above calculation is the 
inclusion of the heavy quark PDF contribution which re- 
sums logs of as ln(/j, 2 /m 2 ). An obvious question is when 
do we need to consider such terms, and how large are 
their contributions? The answer is illustrated in Fig. [4] 
where we compare the DGLAP evolved PDF /q(x,ju) 
with the single splitting perturbative result /q(;t,/x). 

The DGLAP PDF evolution sums a non-perturbative 
infinite tower of logs which are contained in <jt,o while 
the o~suB contribution removes the perturbative single 
splitting component which is already included in the 
o'nlo contribution. Hence, at the PDF level the differ- 
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Figure 5: F% for x — 0.1 for NLO DIS heavy quark pro- 
duction as a function of Q. We display calculations us- 
ing the ACOT, S-ACOT, Fixed-Flavor Number Scheme 
(FFNS), and Zero-Mass Variable Flavor Number Scheme 
(ZM-VFNS). The ACOT and S-ACOT results are virtu- 
ally identical. 



ence between the heavy quark DGLAP evolved PDF /q 
and the single-splitting perturbative /q will indicate the 
contribution of the higher order logs which are resummed 
into the heavy quark PDF. Here, /q = / s <£> Pg^Q rep- 
resents the PDF of a heavy quark Q generated from a 
single perturbative splitting. 

For (j, ~ m we see that /q and /q match quite closely, 
whereas they differ significantly for \i values a few times 
m. While the details will depend on the specific pro- 
cess, in general we find that for /^-scales a few times m 
the terms resummed by the heavy quark PDF can be 
significant. Additionally, the difference between Jq and 
Jq will be reduced at higher orders as more perturbative 
splittings are included in /q. 

Note that these scales are much lower than one might 
estimate using the naive criterion ^ ln(/i 2 /m 2 ) ~ 1; in 
particular, the ACOT calculation often yields reduced in- 
dependence as the quark dominated ctlo contributions 
typically have behavior which is complementary to the 
gluon-initiated ctnlo terms. 



B. S-ACOT 

In a corresponding application, it was observed that 
the heavy quark mass could be set to zero in certain 
pieces of the hard scattering terms without any loss of 
accuracy. This modification of the ACOT scheme goes 
by the name Simplified-ACOT (S-ACOT) and can be 
summarized as follows [8|. 

S-ACOT: For hard-scattering processes 
with incoming heavy quarks or with inter- 
nal on-shell cuts on a heavy quark line, the 
heavy quark mass can be set to zero (in = 0) 
for these pieces. 



If we consider the case of NLO DIS heavy quark pro- 
duction, this means we can set m = for the LO terms 
(QV — > Q) as this involves an incoming heavy quark, 
and we can set m = for the SUB terms as this has an 
on-shell cut on an internal heavy quark line. Hence, the 
only contribution which requires calculation with m re- 
tained is the NLO gV — > QQ process. Figure [5] displays 
a comparison of a calculation using the ACOT scheme 
with all masses retained vs. the S-ACOT scheme; as 
expected, these two results match throughout the full 
kinematic region. 

It is important to note that the S-ACOT scheme is 
not an approximation; this is an exact renormalization 
scheme, extensible to all orders. 



C. ACOT and x-Rescaling 



As we have illustrated in Sec. Ill Al above, in the limit 
Q 2 3> m 2 the mass simply plays the role of a regula- 



tor. In contrast, for Q 



the value of the mass 



is of consequence for the physics. The mass can enter 
dynamically in the hard-scattering matrix element, and 
can enter kinematically in the phase space of the process. 
We will demonstrate that for the processes of interest 
the primary role of the mass is kinematic and not dy- 
namic. It was this idea which was behind the original 
slow-rescaling prescription of [9( which considered DIS 
charm production (e.g., 7c — > c) introducing the shift 



X = x 






(3) 



This prescription accounted for the charm quark mass 
by effectively reducing the phase space for the final state 
by an amount proportional to (m c /Q) ■ 

This idea was extended in the x-scheme by realiz- 
ing that in addition to the observed final-state charm 
quark, if the beam has a charm-flavor quantum number 
of zero (such as a proton beam) then there is also an 
anti-charm quark in the beam fragments because all the 
charm quarks are ultimately produced by gluon splitting 
(g — > cc) into a charm pair|3 For this case the scaling 
variable becomes 



X = x 



2m c 



(4) 



7 If the beam has non-zero charm-flavor quantum number, such 
as a £)-meson, this argument would be incorrect. Technically, 
X-scaling violates factorization as we are presuming the mass of 
the beam fragments; if we perform a thought experiment with 
a beam of D-mesons, charm quark need not be associated with 
an anti-charm quark. 



This rescalingis implemented in the ACOT x scheme, 
for example (lol4ll|. The factor (1 + (2m c ) 2 /Q 2 ) repre- 
sents a kinematic suppression factor which will suppress 
the charm process relative to the lighter quarks. Addi- 
tionally, the x-scaling ensures the threshold kinematics 
(W 2 > 4m 2 + M 2 ) are satisfied; while it is important to 
satisfy this condition for large x, this may prove too re- 
strictive at small x where the HERA data are especially 
precise^ 

To encompass all the above results, we can define a 
general scaling variable x(n) as 



X(n) 



Q 



(5) 



where n = {0,1,2}. Here, n = corresponds to the 
massless result without rescaling, n = 1 corresponds to 
the original Barnett slow-rescaling, and n = 2 corre- 
sponds to the x-rescaling. 



D. Phase Space (Kinematic) & Dynamic Mass 

We now investigate the effects of separately varying 
the mass entering the x( n ) variable taking into account 
the phase space constraints and the mass value entering 
the hard scattering cross section a(m). We call the for- 
mer mass parameter "phase space (kinematic) mass" and 
the latter "dynamic mass'uJ 

In Fig. [Ha|we display F£(x,Q) vs. Q. The family of 
3 curves shows the NLO ACOT calculation with x( n ) 
scaling using a zero dynamic mass for the hard scatter- 
ing. We compare this with Fig. I6bl which shows F£(x, Q) 
in the NLO ACOT scheme using a fixed n — 2 scal- 
ing, but varying the mass used in the hard-scattering 
cross section. The upper (cyan) curves use a non-zero 
dynamic mass [a(m c = 1.3)] and the lower (purple) 
curves have been obtained with a vanishing dynamic 
mass [a(m c = 0)]. We observe that the effect of the 
'dynamic mass' in a(m c ) is only of consequence in the 
limited region Q > m, and even in this region the ef- 
fect is minimal. In contrast, the influence of the phase 
space (kinematic) mass shown in Fig. |Ba] is larger than 
the dynamic mass shown in Fig. [6b] To highlight these 
differences, we scale the curves in Fig. [7] by the massless 
n = 2 scaling result and plot bands that represent the 
variation of the dynamic and kinematic masses. 



8 We sketch the relevant kinematics in Appendix [A3] 

9 Note that the finite mass terms (m 2 /Q 2 ) n in a(m) receive con- 
tributions from both, masses in the heavy quark propagators 
and masses in the phase space. Still we refer to them as dynamic 
mass terms and show that they are numerically less important 
than the mass terms in the slow rescaling variable x(n) which 
are of purely kinematic origin. 



In conclusion, we have shown that (up to 0(as)) the 
phase space mass dependence is generally the dominant 
contribution to the DIS structure functions. Assuming 
that this observation remains true at higher orders, it is 
possible to obtain a good approximation of the structure 
functions in the ACOT scheme at NNLO and N 3 LO us- 
ing the massless Wilson coefficients together with a non- 
zero phase space mass entering via the ^(^-prescription. 



III. ACOT SCHEME BEYOND NLO 

We have shown using the NLO full ACOT scheme 
that the dominant mass effects are those coming from 
the phase space which can be taken into account via a 
generalized slow-rescaling x(n)-prescription. Assuming 
that a similar relation remains true at higher orders, one 
can construct the following approximation to the ACOT 
result up to N 3 LO (0(a|)): 

ACOT[C(a° +1+2+3 )] ~ 

ACOT[0(c4 +1 )] + ZM-VFNS x(ll) [C(a| +3 )] (6) 

In this equation, "ACOT" generically represents any vari- 
ant of the ACOT scheme (ACOT, S-ACOT, S-ACOT x ); 
for the results presented in Sec. IIV1 we will use the fully 
massive ACOT scheme with all masses retained out to 
NLO. The ZM-VFNS x(ll) term uses the massless Wilson 
coefficients at 0(aa 2 s ) and 0(aag) with the specified 
Y(n)-scalingP^I Sample processes which contribute at 
this order are displayed in Fig. [H] 

We use the ZM-VFNS X („) result in Eq. © to ap- 
proximate the higher-order terms because not all the 
necessary massive Wilson coefficients at 0(aa 2 s ) and 
C(aa|) have been computed. There has been a calcu- 
lation of neutral current electroproduction (equal quark 
masses, vector coupling) of heavy quarks at this order 
by Smith & VanNeerven [U] in the FFNS which could 
be used to obtain the massive Wilson coefficients in the 
S-ACOT scheme by applying appropriate collinear sub- 
traction terms. However, for the original ACOT scheme 
it would then still be necessary to compute the massive 
Wilson coefficients for the heavy quark initiated subpro- 
cess at 0(aa 2 s ). See Refs. Il2l . \J4\ for details. 

Using the result of Ref. |13j , Thorne and Roberts de- 
veloped an NLO VFNS 0,111, and an improved NNLO 
formulation was presented in Ref. [ly]- The FONLL 
formalism was outlined in Ref. |17| and this was used 
to construct matched expressions for structure functions 
to NNLO [lji] ; implications of these results in the con- 
text of the NNPDF analysis were presented in Ref. [19| . 



10 In Sec. II. A. 2 we demonstrated that the ACOT calculation re- 
duces to the ZM-VFNS result in the massless limit. We will 
address the choice of the x("-)-rescaling in the Sec. III. A. 
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the scaling fixed n = 2 and compare the effect of varying the dynamic mass in the Wilson coefficient. The upper (cyan) curve 
uses a non-zero dynamic mass [a(m = 1.3)] and the lower (purple) curve uses a zero dynamic mass [a(m = 0)]. 
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Figure 6: Comparison of phase space (kinematic) & dynamic mass effects 
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Figure 7: Comparison of kinematic & dynamic mass effects for F£(x, Q) vs. Q for the NLO ACOT calculation for 

x = {10 -1 , 10 -3 , 10 -5 } (left to right). The curves are scaled by the massless n — 2 result. The wider (yellow) band 

represents the variation of the kinematic mass of Fig. [Sal note this band extends down to a ratio of 1.0. The narrower 

(blue) band is overlaid on the plot and represents the variation of the dynamic mass of Fig. I6bl 



An overview and comparison of these analyses was pre- 
sented in the 2009 Les Houches report [2CJ]. More re- 
cently, an NNLO S-ACOT-^ calculation was developed 
in Refs. fl2i. Il-4l| . For charge current case massive calcu- 
lations are available at order 0(aas) |2ll — [23j | and partial 
results at order O{ao? s ) [2J|. Comparative analyses of 
these schemes are under investigation; however, this is 
beyond the scope of this paper. 

Here, we argue that the massless Wilson coefficients at 
0{aa 2 s ) together with a x(n)-prescription provide a very 
good approximation of the exact result. At worst, the 
maximum error would be of order 0(aa 2 > x [m 2 /Q 2 ]). 
However, based on the arguments of Sec. HID I we expect 
the inclusion of the phase space mass effects to contain 
the dominant higher order contributions so that the ac- 



tual error should be substantially smaller. 

The massless higher order coefficient functions for the 
DIS structure function F 2 via photon exchange can be 
found in Refs. jl!^ for 0(a\), Refs. [H-{30| for 0{a%), 
and Ref. (3l| for 0(a|). For our numerical code we have 
used the x-space parameterization provided in Refs. |32l . 
13 for C(a|), and Refs. (HHl for G(a%). 

The expressions for the structure function Fl have 
been calculated in Refs. HHH] for C(a|), and Ref. HU 
for 0(a 3 s ). In our FORTRAN code we have used the 
x-space parameterization provided in Refs. [32l . |36| for 
C(a|) and Ref. [H| for C(a|). 

In order to calculate the inclusive structure functions 
F 2 and F L in the ZM-VFNS X using these Wilson coeffi- 
cients, plus- and delta-distributions have to be evaluated 
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Figure 8: Sample Feynman diagrams contributing to DIS 
heavy quark production (from left): LO 0(oP s ) quark- 
boson scattering QV — > Q, NLO 0(ag) gluon-boson 
scattering gV — > QQ, NNLO 0{a%) boson-gluon scat- 
tering gV — > gQQ, and N 3 LO C(a|) boson-gluon scat- 
tering gV -> ggQQ. 



which is in principle straightforward. However, for the 
implementation of the slow-rescaling prescription it is 
necessary to decompose the Wilson coefficients into the 
contributions from different parton flavors. This step is 
non-trivial at 0(a 2 s ) and beyond, and we therefore pro- 
vide some details of our calculation in the Appendix [Bj 



A. Choice of x(^)-R.escaling 



IV. RESULTS 

We now present the results of our calculation extend- 
ing the ACOT scheme to NNLO and N 3 LO. As out- 
lined in Eq. ((6]), we will use the fully massive ACOT 
scheme for the LO and NLO contributions, and com- 
bine this with the ZM-VFNS supplemented with the %- 
rescaling prescription to approximate the higher order 
terms. We will use the QCDNUM program [H] with 
the VFNS evolved with the DGLAP kernels at NNLO 
to generate our PDFs from an initial distribution based 
on the Les Houches benchmark set [39| ; this ensures that 
our heavy quark PDFs are consistently evolved so that 
the heavy quark initiated LO terms properly match the 
corresponding SUB contribution. At NNLO the proper 
matching conditions across flavor thresholds introduces 
discontinuities in the PDFs which are incorporated in 
the QCDNUM program; we discuss this in detail in Ap- 
pendix C. We choose m c = 1.3 GeV, m& = 4.5 GeV, 
as(M z ) = 0.118. We note that the QCDNUM ZM- 
STFN package has the massless Wilson coefficients com- 
puted up to N 3 LO; we cross checked our implementation 
of ACOT in the massless limit with QCDNUM, and they 
agree precisely. 



We now consider our choice for the appropriate gen- 
eralized x(n)-rescaling variable. 

In Table U we display the various rescalings of £ for 
the LO 7Q —> Q process and the NLO 73 — > QQ pro- 
cess. The "general" result is obtained by working out the 
detailed kinematics for the corresponding process |37| . 

The factor r\ is the rescaling due to the hadronic mass 
M; notice that this factors out from the partonic mass 
dependence as it should [lj. For details see Appendix lAl 

The LO case with full massive kinematics has been 
computed in Ref. [37|. In the limit where the initial 
mass is small (rti\ — > 0), we recover the Barnett J9fl slow- 
rescaling result. Additionally, we obtain the curious re- 
sult that for a neutral current equal mass case (mi = 1112) 
the rescaling is this same factor. 

For the NLO gluon-induced process, the interpreta- 
tion of the rescaling is straightforward; the phase space 
is simply suppressed by the total invariant mass of the 
final state (m\ + 111,2) compared to the scale Q. For 
the charged current case where we neglect mi, we again 
obtain the standard rescaling factor. However, for the 
neutral current case (mi = 777,2) we obtain a rescaling 
factor which is analogous to the ^-scaling factor. 

For the purposes of this study, we will vary the phase 
space mass using the x( n ) rescaling with n = {0, 1,2}. 
While 77 = corresponds to the massless case (no rescal- 
ing) , it is not obvious whether 77 = 1 or 77 = 2 is the pre- 
ferred rescaling choice for higher orders. Thus, we will 
use the range between 77 = 1 and 77 = 2 as a measure of 
our theoretical uncertainty arising from this ambiguity. 



A. Effect of x(7i)-Scaling 

In Figures [9a] and I9bl we display the structure func- 
tions F2 and Fl, respectively, for selected x values as 
a function of Q. Each plot has three curves which are 
computed using 77-scalings of {0,1,2}. We observe that 
the effect of the 71-scaling is negligible except for very 
small Q values. This result is in part because the heavy 
quarks are only a fraction of the total structure function, 
and the effects of the n-scaling are reduced at larger Q 
values. 

In Fig. [10] we magnify the small Q region of Fl of 
Fig. [9b] for x = 10~ 5 , where the effects of using different 
scalings are largest. We can see that for inclusive observ- 
ables, the 77 = 1 and 77 = 2 scalings give nearly identical 
results, but they differ from the massless case (77 = 0). 
This result, together with the observation that at NLO 
kinematic mass effects are dominant, suggests that the 
error we have in our approach is relatively small and 
approximated by the band between n = 1 and 77 = 2 
results. 



B. Flavor Decomposition of xi n ) Scaling 

We can investigate the effects of the x(n)-scaling in 
more details by examining the flavor decomposition of 
the structure functions. 

In Figures fllal and rilbl we display the fractional contri- 
butions of quark flavors to the structure functions F^l 



General 



Q 1 -mf+mf 2 +A[-Q A ,mj,m%] 
2Q^ 



( milmj 



mi = 



V 1 + 



mi = rri2 = m 



1 + 



•J 1 + ^ 



X-scheme: 



'/ 



1 -i- < 2m > 



»J 1 + ^ 



Table I: The massive rescaling factor for the LO quark-initiated process (Vqi —> (72), and the NLO gluon-initiated 
process (Vg — > gi^)- The quarks qi^ have mass mi. 2, respectively, and T^ represents the vector boson; j/Z 
for neutral current processes (mi = TO2), and W for charged current processes (mi 7^ Ttiz). r\ is the scaling 
factor which depends on the hadronic mass M; see Appendix [3] for details. The triangle-function is defined as: 

A[a,6,c] = y/a 2 + b 2 + c 2 - 2(ab + bc + ca). 
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Figure 9: F 2 .l vs. Q at N 3 LO for fixed x = {lO -1 , 10~ 3 , 10~ 5 } (left to right). The three lines show the scaling 
variable: n — {0, 1, 2} (red, green, blue). We observe the effect of the n-scaling is negligible except for very small Q 

values. 
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Figure 10: Enlargement of Fig. [9b] for x = 1CP 5 showing 
the small Q region. Here we can distinguish plots for dif- 
ferent scalings; from top to bottom we have n = {0, 1, 2} 
(red, green, blue). 



fined in appendix [B] in Eqs. ()B"T1) and (|B2|) . We ob- 
serve the n-scaling reduces the relative contributions of 
charm and bottom at low Q scales. For example, with- 
out any n-scaling (n = 0) we find the charm and bot- 
tom quarks contribute an unusually large fraction at very 
low scales (Q ~ m c ) as they are (incorrectly) treated as 
massless partons in this region. The result of the dif- 
ferent n-scalings (n = 1,2) is to introduce a kinematic 
penalty which properly suppresses the contribution of 
these heavy quarks in the low Q region. In the following, 
we will generally use the n = 2 scaling for our compar- 
isons. 



C. F2,l Initial-State Flavor Decomposition 



for selected n-scaling values as a function of Q. Fla- 
vor decomposition of inclusive structure functions is de- 



In Figures I12al and I12bl we display the fractional con- 
tributions for the initial-state quarks (i) to the structure 
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Figure 11: Effect of x("-)-scaling for n = {0,1,2} (left to right) at N 3 LO for fixed x = {10 -3 }. Reading from the 
bottom we have fractional contribution for each quark flavor to F^ l JFi.l vs. Q from {w, <i, s, c, b} (green, blue, cyan, 

magenta, pink). 



functions F% and .FlCJ respectively, for selected x val- 
ues as a function of Q; here we have used n = 2 scal- 
ing. Reading from the bottom, we have the cumulative 
contributions from the {g, u, d, s, c, b}. Although this de- 
composition is not physically observable, it is instructive 
to see which PDFs are dominantly influencing the result. 
We observe that for large x and low Q the heavy flavor 
contributions are minimal. For example, for x = 10 _1 
we see the contribution of the u-quark comprises ~ 80% 
of the i<2 structure function at low Q. In contrast, at 
x = 10 -5 and large Q we see the F2 contributions of 
the w-quark and c-quark are comparable (as they both 
couple with a factor 4/9), and the d-quark and s-quark 
are comparable (as they both couple with a factor 1/9). 
It is notable that the gluon contribution to Fl is sig- 
nificant. For x = 10" 1 this is roughly 40% throughout 
the Q range, and can be even larger for smaller x values. 



D. F2,l Final-State Flavor Decomposition 



In Figures I13al and I13bl we display the fractional con- 
tributions for the final-state quarks (J) to the structure 
functions F2 and Fl , respectively, for selected x values as 
a function of Q; here we have used n = 2 scaling. Read- 



ing from the bottom, we have the cumulative contribu- 
tions from the {u, d, s, c, b}. Again, we observe that for 
large x and low Q the heavy flavor contributions are min- 
imal, but these can grow quickly as we move to smaller 
x and larger Q. 



E. Comparison of LO, NLO, NNLO, N 3 LO 



Fractional decomposition of "initial-state" structure functions is 



understood as Fl 



= 1 r 2,L- 



In Figure I14al we display the results for F2 vs. Q 
computed at various orders. For large x (c.f. x — 0.1) we 
find the perturbative calculation is particularly stable; 
we see that the LO result is within 20% of the others at 
small Q, and within 5% at large Q. The NLO is within 
2% at small Q, and indistinguishable from the NNLO 
and N 3 LO for Q values above - 10 GeV. The NNLO 
and N 3 LO results are essentially identical throughout the 
kinematic range. For smaller x values (10~ 3 , 10 -5 ) the 
contribution of the higher order terms increases. Here, 
the NNLO and N 3 LO coincide for Q values above ~ 
5 GeV, but the NLO result can differ by - 5%. 

In Figure [Mb] we display the results for Fl vs. Q com- 
puted at various orders. In contrast to F%, we find the 
NLO corrections are large for Fl ; this is because the LO 
Fl contribution (which violates the Callan-Gross rela- 
tion) is suppressed by (to 2 /Q 2 ) compared to the dom- 
inant gluon contributions which enter at NLO. Conse- 
quently, we observe (as expected) that the LO result 
for Fl receives large contributions from the higher or- 
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Figure 12: Fractional flavor decomposition of "initial-state" F\ l JFi,l vs. Q at N 3 LO for x = {10" 1 , 10 -3 , 10~ 5 } 
(left to right) for n = 2 scaling. Reading from the bottom, we plot the cumulative contributions to F 2 ,l from 

{g, u, d, s, c, 6}, (red, green, blue, cyan, magenta, pink). 



der termso Essentially, the NLO is the first non-trivial 
order for Fl, and the subsequent contributions then 
converge. For example, at large x (c.f. x = 0.1) for 
Q ~ 10 GeV we find the NLO result yields - 60 to 
80% of the total, the NNLO is a ~ 20% correction, and 
the N 3 LO is a ~ 10% correction. For lower x values 
(10 -3 , 10 -5 ) the convergence of the perturbative series 
improves, and the NLO results is within ~ 10% of the 
N 3 LO result. Curiously, for x = 10~ 5 the NNLO and 
N 3 LO roughly compensate each other so that the NLO 
and the N 3 LO match quite closely for Q > 2 GeV. 

While the calculation of Fl is certainly more chal- 
lenging, examining Fig. [T] we see that for most of the 
relevant kinematic range probed by HERA the theoret- 
ical calculation is quite stable. For example, in the high 
Q 2 region where HERA is probing intermediate x values 
(x ~ 10~ 3 ) the spread of the x( n ) scalings is small. The 
challenge arises in the low Q region (Q ~ 2 GeV) where 
the x values are ~ 10~ 4 ; in this region, there is some 
spread between the various curves at the lowest x value 
(~ 10~ 5 ), but for x ~ 10 -3 this is greatly reduced. 



12 Because we use the fully massive ACOT scheme to LO and NLO, 
the LO result in Fig. ll4bl contains the (m 2 /Q 2 ) helicity- violating 
contributions ~ O(a^); hence, it is non-zero. In the S-ACOT 
scheme, the LO result for Fl vanishes, but the NLO result is 
comparable to the NLO ACOT result. 



CONCLUSIONS 



We extended the ACOT calculation for DIS struc- 
ture functions to N 3 LO by combining the exact ACOT 
scheme at NLO with a x(n)-rescaling; this allows us 
to include the leading mass dependence at NNLO and 
N 3 LO. Using the full ACOT calculation at NLO, we 
demonstrated that the heavy quarks mass dependence 
for the DIS structure functions is dominated by the kine- 
matic mass contributions, and this can be implemented 
via a generalized x(n)-rescaling prescription. 

We studied the F 2 and Fl structure functions as a 
function of x and Q. We examined the flavor decom- 
position of these structure functions, and verified that 
the heavy quarks were appropriately suppressed in the 
low Q region. We found the results for F 2 were very 
stable across the full kinematic range for {x, Q}, and 
the contributions from the NNLO and N 3 LO terms were 
small. For Fl, the higher order terms gave a propor- 
tionally larger contribution (due to the suppression of 
the LO term from the Callan-Gross relation): neverthe- 
less, the contributions from the NNLO and N 3 LO terms 
were generally small in the region probed by HERA. 

The result of this calculation was to obtain precise pre- 
dictions for the inclusive F 2 and Fl structure functions 
which can be used to analyze the HERA data. 
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Figure 13: Fractional contribution for each quark flavor to F| l /F2 ! l vs. Q at N 3 LO for fixed x = {10 1 , 10 3 , 10 5 } 
(left to right). Results are displayed for n = 2 scaling. Reading from the bottom, we have the cumulative contributions 

from the {w, d, s, c, b} (green, blue, cyan, magenta, pink). 
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to right) for n = 2 scaling. 



Appendix A: Kinematic Relations 

1. Target Mass Contributions 

In the DIS process, the effect of the target mass (M) 
on the scaling variable is a multiplicative correction fac- 



tor 



2.7- 



-, , 4x 2 M 2 
1 H 77T- 



AI^Q 



fxM\ 



(Al) 
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Figure 15: O(a s ) - 7*% -¥ q t . Contributes to C™ q (and Figure 16: 0{a 1 s ) - 7*5, -> ^.g. Contributes to C£ s q (and 

hence to C* „) but not to CT S „. This contribution does 



hence to C a , 



but not to C** q . 



but not to C%* q . 

not depend on rif. 



This is used in Table U to modify the scaling variable [J, 



2. Barnett Scaling 

If we consider the charged-current DIS process for 
charm production, this takes place via the subprocess 
W + (q)s(t;P) — > c(k). If we impose 4-momentum con- 



servation, we have (q + t;PY 



q 2 = -Q 2 and x = Q 2 /(2r 
"slow rescaling" relation [£ 



Defining 



• g) , we obtain the traditional 



x 1 



Q 2 



which was used in Eq. ([3]). 



W constraints 



If we compute the invariant mass W of a boson of mo- 
mentum g scattering from a light parton a of momentum 
p a = £P, we find |li| 



W = (p« + g) 2 = <9 2 (£/z - l) 



(A2) 



If the partonic final state has a minimum invariant mass 
Wmin = 4m 2 , then £ is constrained by 

1 > £ > X > x (A3) 

where \ = a?(l + 4m 2 /Q 2 ). This is the relation used 
in Eq. IQ}. This choice will ensure W > W m i n is satis- 
fied. While this constraint is important in the large x 
region, this may be too restrictive in the small x region — 
especially as this is the region where the HERA data is 
very precise. 



Appendix B: Decomposition of the Wilson 
coefficients 



In this appendix we present the decomposition of the 
Wilson coefficients used to implement the scheme. We 
will need to decompose the structure function F in terms 
of the individual partonic contributions, 



F = J2J2 Fl 



(Bl) 




Figure 17: 0{a\\ 
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qjQj- 



where the indices i and j represent initial and final-state 
partons respectively (see captions of Figs. rT5HlM)) . More 
specifically, i = denotes a gluon and i,j = 1, 2, 3, . . . 
denotes u, d, s, ... quarks and anti-quarks. A top quark 
PDF (i = 6) is not included in this study. 

Let us consider the heavy quark structure functions 
F£ L as an example. This is obtained by requiring that 
there is a charm in the initial state while summing 
over the final-state flavors up to and including charm 
in Eq. (|B1[) , or by requiring that there is a charm in the 
final-state and summing over the initial flavors up to and 
including charm. Thus, we obtain: 



pc = J2 pli + H pAj 



F 4 



(B2) 



The case where the initial and final-state are both charm 
quarks (F 44 ) has been written explicitly in the equation 
to avoid double counting this contribution £j The first 
sum in Eq. (|B2|) includes cases, as in Fig. [201 where the 
incoming quark is a light quark while the charm quark 
is one of the quarks in the quark anti-quark pair. 

In order to obtain the required decomposition, there 
are some manipulations that need to be performed 
to transform from the singlet (s), non-singlet (ns), 
and purely-singlet (ps) structure function combinations 
found in the literature into individual partonic compo- 
nents. 

The general expression for the structure function is 
given by: 



l F 

J- n 



Ins ® C% q + (e 2 ) (g s ® C s a<q + g® C a . g ) 



(B3) 



i=0 j=l 



13 Note that in our decomposition, diagrams with a bottom quark 
in the initial or final state, contribute to the bottom structure 
function, even in the presence of a charm quark. 
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Figure 18: 0(a 2 s ) - 
(and hence to C s a ) 



~f*qi — > qtgg. Contributes to C™ 
but not to C*p s . This part is inde- 
pendent of rif . 



Figure 19: C(a|) - j*g -> q 3 q 3 g. 
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(b) Contribution propor- 
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Figure 20: 0(a%) - 7*& -* qiq 3 q 3 . 




Figure 21: (9(a|) - 7*17, 



ft; 953 ■ 



Contribution to C, 
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Figure 22: 0(a%) - 7*3 -> <7j9j33- 

where a = {2, L}, and 

"/ 
fe8 = £(e? - {e 2 ))qf 

8=1 

<7s = JI^' < li = ^ + Qi 

^£ e2 - 






n/ 



(B4) 



and C" s „, Ci! „, C a are the Wilson coefficients. From 

Q,t/ ~ liy' u, ;y 

Eq. (|F34[) one can extract the contribution from a single 
initial-state quark as: 



l F 



4 ® [e? C£ 



>C£TJ 



where C% s q is 



C ps 

a,q 



c, 



a . q 



c, 



a .</ 



(B5) 



(B6) 



To further decompose Eq. (|B5|) into the different final- 
state contributions, we examine the diagrams that con- 
tribute to the non-singlet and purely-singlet coefficients. 
Diagrams in which the photon couples to the incoming 
quark contribute to C™ q (Figs. El EEBJ EDU HE etc.), 
whereas the diagrams where the photon does not couple 
to the incoming quark contribute to C^ s q ; these contribu- 
tions appear for the first time at <D(a 2 s ) in Figs. l20all23al 
Separating out the final-state quark from Eq. (|B5[) we 
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(a) Contribution propor- 
tional to rif for Ca^q- 



Eq. (|B7[) since there are only purely-singlet contri- 
butions: 



x- l FV=g®Ue 2 )^C atg (j) 
-{e 2 )U-VC a , g (j-l) 



(B8) 




(b) Contribution propor- 
tional to nj for C™ . 



Figure 23: C(a|) - 7*ft -> qiq^g- 



Tnnr^^f 
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Figure 24: 0(a|) - 7*3 -» qjq 3 q k qk- 



obtain: 



3.-1^ - ~+ < 



q[ ® < e 



C(n/ = 0) *, 



+ qj)-Qj-i) (B7) 



+ {e 2 ) {j) cz: q {j)-{e 2 ) {j - x) cz: q {j--t) 



We have introduced <5jj in the non-singlet contribution 
to account for contributions in which the photon couples 
to the initial and final-state quark. When this is not the 
case, (i.e., in all purely-singlet contributions and in non- 
singlet contributions such as the ones in Fig. I20b| , the 
difference of the coefficient functions with rif = j and 
rif = j — 1 flavors is taken. 
Some comments are in order: 

• We have verified analytically and numerically that 
one recovers Eq. (|B5|) when summing over the final 
state quark partons (j = 1, . . . ,n/) in Eq. (|B7|) . 

• The corresponding decomposition for the gluon- 
initiated subprocesses is simpler than the one in 



• We remark that the decomposition in Eq. (|B7[) also 
includes the contributions from virtual diagrams to 
the Wilson coefficients. As has been discussed in 
the literature [40J , such a decomposition is ambigu- 
ous at 0(a 2 s ) and beyond due to the treatment of 
heavy quark loops contributing to the light quark 
structure functions. However, numerically the am- 
biguous terms are small and it is standard to ana- 
lyze the heavy quark structure functions F£ L and 
F% L in addition to the inclusive structure functions 
F2.L without any further prescription. 

For the general neutral current case (including Z- 
boson exchange), the electromagnetic couplings should 
be replaced by electroweak couplings as follows: 



■ 2 ei v e v q xz + {vl + al){v 2 q + a\) X \ (B9) 



where 



Vf 



Tf 



2Q f sin 2 6>v 



'/ 



— J 3 



(BIO) 



are the standard (axial-)vector couplings of the Z-boson 
to the leptons (/ = e) and quarks (/ = q). Furthermore, 
Xz is the ratio of the Z-boson propagator with respect 
to the photon propagator including additional coupling 
factors: 



Xz 



G F M 2 Z 



Ml 



2y/2iza e m Q 2 
Finally, the average squared charge is modified as 



(Bll) 



(e 2 ) 



2\("/) 



1 _. 

a+ ("/) = — E< 



(B12) 



Appendix C: Matching Across Heavy Flavor 
Thresholds 

As we compute at higher orders, we find the matching 
conditions of the PDFs become discontinuous at 0(a 2 s ) 
(NNLO), and the matching of the MS a s (n) becomes 
discontinuous at 0{a 3 s ) (N 3 LO). 

While the discontinuities in the PDFs and a s (which 
are unphysical quantities) persist at all orders, physical 
observables (such as cross sections and structure func- 
tions) will match across thresholds up to the computed 
order of the perturbation theory; for example, a physical 
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observable in an N- flavor and an (N + l)-flavor scheme 
will match up to higher order terms when computed to 
order a^ in the perturbation expansion: 



r N 



r N+l 



+ C(af +i ) 



As it is not immediately obvious how the discontinu- 
ities cancel order- by-order, we shall examine a NNLO 
numeric case, and also a simple analytic example. 



1. Discontinuities across the flavor transition. 

To illustrate the behavior of the discontinuities, we 
will work at NNLO where the DGLAP evolution and 
the flavor-threshold boundary conditions have been com- 
puted and implemented O Since \x = m c is often used for 
the initial evolution scale, we will focus on the transition 
from Np = 4 to Np = 5 flavors at /j, = mj,. 

The matching conditions across flavor thresholds can 
be summarized as |4ll 



f 

J a 



N+l 



A n 



fx 



N 



(CI) 



where f N and f N+1 are the PDFs for N and N + 1 fla- 
vors, and A a b can be expanded perturbatively. In the 
VFNS for n < 77ib, the 6-quark PDF is zero and the 
gluon PDF is finite and positive. Using Eq. (|C1|) for 
\x > ?77f,, we find the 6-quark is negative for fj, ~ mj,, 
and it becomes more negative as we move to smaller x. 
In contrast, the gluon has a positive discontinuity as it 
must to ensure the momentum sum rule is satisfied. 



2. The 6-quark flavor transition. 

Although these discontinuities are too small to be no- 
ticeable in the figures of Sec. IV, in Figure [5S] we have 
magnified the axes so the discontinuities are visible. 
Here, we display F2 and Fl for a selection of x- values. 

The first general feature we notice in Fig. [25] is that 
the size of the discontinuity generally grows as we go to 
smaller x values. This is consistent with the fact that 
the discontinuity computed by Eq. (|C1|) also grows for 
smaller x. We display the results for a selection of 77- 
scaling values; note that the uncertainty arising from the 
discontinuity is typically on the order of the difference 
due to the choice of scaling. 

Another feature that is most evident for the series of 
Fl plots (Fig. I25b[) is that the discontinuity can change 
sign for different x values. This can happen because the 
mix of quark and gluon initiated terms is changing as a 
function of x. 



This observation is key to understanding how the (un- 
physical) PDFs may have a relatively large discontinu- 
ity, while the effect on the physical quantities (such as 
a and F2.L) is moderated. Because physical quantities 
will contain a sum of gluon and quark initiated contri- 
butions, and because the discontinuity of the quark and 
gluon PDFs have opposite signs, the discontinuities of 
the quark and gluon PDFs can partially cancel so that 
the physical quantity may have a reduced discontinuity. 

This discontinuity, in part, reflects the theoretical un- 
certainty of the perturbation theory at a given order. 
As we compute the physical observables to higher and 
higher orders, this discontinuity will be reduced even 
though the discontinuity in the PDFs and a s remain. 
We will demonstrate this mechanism in the following. 



3. A "Toy" Example at NLO" 

We now illustrate how the cancellation of the quark 
and gluon PDF discontinuities work analytically using a 
"toy" calculation. 

Expanding Eq. (|C1|) in the region of /j = 777;, we have: 



f! = { 



= 



^(L 



J a 



2- 



P aa {L 



,) + 0{a 2 s ))®f l 
7 ) + 0(a 2 s )}®f ! 



where L 
terms a 



= ln{[i 2 /ml). 



ij 



(C2) 

It happens that the constant 
in Eq. ([C2|) are zero at 0{a\) in the WS 
scheme; this is not due to any underlying symmetry, and 
in fact at 0{a 2 s ) these terms are non-zero. Because a^ 
are zero, if we perform the matching at \x = m^, we find 
that the gluon PDF is continuous f^(x, rtib) = fg(x, rrib), 
and the bottom PDF starts from zero f^(x,mb) = 0. 



a. If dij was non-zero at 0{a\) 

To illustrate how the discontinuities cancel in the 
ACOT renormalization scheme, we will suppose (for this 
"toy" calculation) that the constant terms (a^) in the 
matching conditions are non- vanishing at order a\; thus, 
the gluon and bottom PDFs will now have 0{a\) dis- 
continuities, but the physical observables computed with 
different Np values will still match up to 0(a 2 s ). 

In the ACOT scheme, the total cross section can be 
decomposed as: o T ot = ^lo+^nlo-^suBj where a LO 
represents 76 — > 6, o^lo represents jg — > bb, and o~sub 
represents the (g —>■ b) ® (76 — > b) "subtraction" contri- 
bution!^ We will now perturbatively compute o~tot in 
the region /i ~ 777^ for both Np = 4 and Np = 5. 



14 At present, the full set of matching conditions and DGLAP ker- 
nels have not been computed at N 3 LO. 



15 Note, we will focus on the gluon-initiated terms, but the demon- 
stration for the quark-initiated pieces is analogous. 
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(b) Fl vs. Q for x = {10 1 , 10 3 , 10 5 } (left to right) for different n scalings. 

Figure 25: Discontinuity for F2, Fp at NNLO in the region of the bottom mass, m& = 4.5 GeV. 



ACOTforN F = 4 



d. Comparison of N F — 5 and Np 



For [a < nib, we have Nf = 4 and ft, = 0; thus, apo 
and (Jsub vanish, and we have: 



N F =A 
'TOT 



ONLO = C l ®.f A Q +0{a 2 s ) 



where C 1 represents the 0(a l s ) process jg — > bb. 



c. ACOTforN F =5 

For /j, > mi,, we have Np = 5 and /& ^ 0. For the 
contributions we have: 



olo 

O'NLO 



C° f b c C° ® {0 + ^P qg (L + a qg )} ® f, 
C 1 ® f g ~ C 1 ® {l + ^-P gg (L + a gg )} I 



&SUB = C ® fg^ q 

a 



/, 



{ 



1 



2tt 



P gg (L 



C°®{^-P qg (L 
* M )}® ft 



% qg)} 



Keeping terms to 0{a\) we have: 



_N F =5 



C l ®ft + 0{al) 



"TOT — a LO + O'NLO ~ ° ' SI) B — - o Jg 

Notice that the discontinuity introduced by a qg in the 
PDFs is canceled by a qg from the SUB contribution Pi 



16 The explicit form of the ACOT subtraction is defined in Sec. IV. C 
(c/., Eq. (36)) of Ref. Q- For an example of the cancellation 
between a xo an d ctsjjb in a more general context, see Ref. |42|| . 



Comparing the Np = 5 and Np = 4 results, we find 



N F =5 
°TOT 



N F =4 
'TOT 



+ 0{a 2 s ) 



so that the total physical results match up the order of 
the perturbation theory. 

In the above illustration, we have retained the log 
terms (L); the cancellation of the logs is ensured in a 
well defined renormalization scheme, and the a^ con- 
stant terms get carried along with the logs and will thus 
cancel order by order. 

Therefore, the discontinuity of the physical quantities 
(o^ ^2,l) reflects the perturbative uncertainty, and this 
will be systematically reduced at higher orders. 
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